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ABSTRACT 

We investigated the probability that an inelastic collision of planetesimals within the Hill 
sphere of the Jovian planets could explain the presence and orbits of observed irregular satel- 
lites. Capture of satellites via this mechanism is highly dependent on not only the mass of 
the protoplanetary disk, but also the shape of the planetesimal size distribution. We per- 
formed 2000 simulations for integrated time intervals ~ 2 Myr and found that, given the 
currently accepted value for the mi nimum mass solar nebula and planetesim al number den- 
sity based upon the lNesvornv et al] d2003l) and lCharnoz & M orbidelli (2003) size distribution 
dN ~ D~ 3 5 dD, the collision rates for the different Jovian planets range between ~ 0.6 
and > 170 Myr -1 for objects with radii, 1 km < r < 10 km. Additionally, we found that 
the probability that these collisions remove enough orbital energy to yield a bound orbit was 
< 10 -5 and had very little dependence on the relative size of the planetesimals. Of these col- 
lisions, the collision energy between two objects was > 10 3 times the gravitational binding 
energy for objects with radii ~ 100 km. We find that, capturing irregular satellites via col- 
lisions between unbound objects can only account for ~ 0.1% of the observed population, 
hence can this not be the sole method of producing irregular satellites. 
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1 INTRODUCTION 

The irregular satellites of the Jovian planets have been the sub- 
ject of much debate since their discovery in the early 1900's. 
Since then there have been approximately 60 found around Jupiter, 
about 60 around Saturn, 20 around Uranus and 10 around Neptun e 
dGladman et all 1200 ll: iHolman et alJl20ol IShe ppard et all 120061) . 
IChamberlin & Yeomansf i 2010h present a current and detailed com- 
pilation of the irregular satellites (e.g. size, orbital characteristics, 
etc.) for the gas giant planets as well as a list of references contain- 
ing these data. 

Although a clear definition of an irregular s atellite remains 
to be determined, we use the same definition as iNesvornv et al] 
d2003h . moons that are far enough from their parent planet that 
the precession o f their orbital plane is primarily controlled by the 
Sun. As stated in Nesv ornv et al] d2003l) . this definition of an irreg- 
ular satellite excludes Neptune's moon Triton. The origin of these 
irregular satellites has yet to be clearly understood and can give 
an insight to the formation of our Solar system. Traditional theo- 
ries of planet formation and evolution fail to produce the observed 
large orbital inclination observed for many of these satellites, some 
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of which even follo w retrograde orbits with r espect to the orbital 
plane of the planets. Ijewitt & Sheppard (20051) review the possible 
modes that possible modes of capturing irregular satellites by the 
giant planets are; 

• Gas Drag: Similar to the Sun and Solar nebula, the gas gi- 
ants are believed to be the result of the gravitational collapse of 
a planetary gas nebula. Therefore, it would be possible for irreg- 
ular satellites to have enough energy dissipated by the planetary 
nebula via gas dra g to remain gravitationally bound to the planet. 
iPollack et al] (1979J) model these phenomena as a possible explana- 
tion for the capture and fragmentation of Jupiter's prograde and ret- 
rograde irregular satellites, but don't provide a clear explanation for 
the irregular satellites orbiting Uranus and Neptune. Furthermore, 
with out the dissipation of planetary nebula, once c aptured, the cap- 
tured object would spiral into the planet ~ 10 yrs dNesvornv et al] 
120031) . 

• Pull Down: Pull down is the sudden mass-growth of the par- 
ent p lanet thereby capturing the satellite (Hepp enheimer & Pored 
1 1977b . This mode of capture occurs when an object enters the Hill 
sphere (region around a planet with radius, Rh = a(m p / Mq) 1 ' 3 , 
where the Sun's gravity is negligible compared the planet's) 
through a Lagrange point and while at this semi-stable orbital point 
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(~ 100 yrs), the mass of the parent planet increases enough for the 
object to become gravitationally bound to the planet. 

• Collisions: Gravitational three body interactions (both collid- 
ing and non-colliding) occurring within the Hill sphere of the planet 
could lead to the capture of one of the objects as well as fragments 
resulting from high en e rgy collisions dColombo & Franklin|[l97ll : 
IWeidenscMllindl2002h . lNesvornv et al.l J2003h calculate the colli- 
sion rate between the existing irregular satellites of Jupiter to be 
< 1 every 4.5 Gyr, implying that their formation hardens back to a 
time when the Solar system was much denser. 

The purpose of this paper is to expand upon the last 
scenario, which has achieved some popularity in recent years. 
iNesvornv et al.l d2004l) that showed that the spectral characteristics 
of the irregular satellite groups around Jupiter and the other plan- 
ets were surprisingly similar. This has been taken to suggest that 
they are all drawn from the same underlying progenitor popula- 
tion, and captured during the early dynamical evolution of the solar 
system. The most detailed scenario for such an evolution is the so- 
called Nice model for the solar system origin dNesvornv et alj2007l ; 
iBottke et al.l20103) . The initial conditions in this model are chosen 
to meet specific benchmarks of a larger scenario involving the ex- 
planation of various dynamical structures in the Solar system, and 
do seem to provide a path to capturing irregular satellites around 
some of the giant planets. Nevertheless, it is also useful to isolate 
the conditions necessary for a specific observable. Thus, in this pa- 
per we wish to evaluate the requirements of a simpler model, in 
which we determine the conditions that allow the irregular satellite 
populations to result from collisional interactions within the tidal 
field of a giant planet. This will also allow us to estimate the size 
of nebula that would match conditions ne cessary to the in situ for - 
mation of Uranus & Neptune proposed by Goldre ich et al. (2004b). 
If irregular satellites could be the remnants of collisions that were 
captured by the planet as a result of said collision, collision rates 
and probabilities of remnants being captured by the planets could 
be used to set limits on the minimum mass Solar nebula (MMSN) as 
well as sizes and number densities of planetesimal seeds that even- 
tually formed the inner planets, asteroid belt, Edgeworth-Kuiper 
belt objects and the Oort cloud. 

The following section, § [2] briefly describes the epoch in the 
formation of the Solar system when the irregular satellites were 
most likely to have been captured. §[3]discusses why it is important 
for us to understand how these satellites were captured by their host 
planets. The presentation of how the data were generated, the moti- 
vation behind it and justification for it are contained in § 14.21 which 
is followed by a detailed description of how various probabilities 
were determined leading up to the calculation of the probability 
that a mass remains bound to a planet, § 15.11 The paper concludes 
with the results of these analyses, §|6] conclusions that can be drawn 
from these results, §|7] and suggested future work, §[8] 



dKokubo & Idal f 1998). During this epoch, the planets in the outer 
Solar system had essentially finished forming, while the terrestrial 
pl anets were still acreting material. A'-Body simulations presented 
in lGoldreich et alj d2004al) and references therein show that for this 
to occur, the Solar nebula must be approximately equal in mass to 
the minimum mass in the inner Solar system and about six times 
the minimum mass of the Solar nebula in the outer solar system. 
This distribution of material is required to ensure the formation of 
Uranus and Neptune in their current location within the age of our 
Solar system (the actual surface densities used were a = 7g cm~ 2 
at 1 AU and a = 1.5 g cm ~ 2 at 25 AU) dGoldreich et al.l2004bl) . 

The epoch directly following oligarchy and thought to be 
the longest period in the evolution of our Solar system is "clean- 
up". At the beginning of this phase of the evolution of our So- 
lar system Jupiter and Saturn are likely to have been completely 
formed, Uranus and Neptune almost completely so. It is likely that 
> 1000 of small bodies were distributed amongst annulii with 
gravitational unstable orbits, to account for the expected distribu- 
tion of ~ 100 of small bodies ejected by Uranus and Neptune 
and a current estimate o f 1 — 10 for the mass of the Oort cloud 
dGoldreich et al.l2004bl) . 



3 IRREGULAR SATELLITES : A WINDOW TO THE 
PAST 

As stated in the first § Q] the irregular moons of the giant planets, 
like Edgeworth-Kuiper belt objects, are thought to be remnants of 
the early Solar system. With this in mind, an accurate description 
of their "capture" by their host planet and possibly the evolution of 
their bound orbits will provide insight to conditions in early Solar 
system. 

As stated in § Q] it is difficult to produce a model where 
the irregular satellites form from the same circumplanetary disk 
from which the regular satellites formed, implying that the irregu- 
lar satellites were captured by the planet through some other means. 
For these satellites to be captured through via collisions, the objects 
collide within the Hill sphere of the planet and this collision must 
dissipate enough energy to result in at least one of objects remain- 
ing gravitationally bound to the host planet. Therefore, we need to 
determine the probability that a planetesimal is located around the 
planet, then a distribution of where around the planet this object is 
most likely to be encountered by a second object. The probability 
of two objects colliding can be calculated from the product of these 
two distributions and a ratio of cross sectional areas (§[6). Finally, 
a collision rate can be calculated by approximating the number of 
planetesimal present near the planets from accepted values of the 
minimum mass Solar nebula and assuming the density of these ob- 
jects is comparable to that of the Earth. 



2 THE EARLY SOLAR SYSTEM: SEEDS FOR 
IRREGULAR SATELLITES 

iGoldreich et"aH d2004bh and references therein give a detailed de- 
scription of the early stages of planetary formation of our Solar 
system. The era most important to this investigation is that of oli- 
garchy since this represents the time just prior to the formation of 
the inner terrestrial planets. At this point, the cores and the gas en- 
velopes of the gas giant planets have already formed. Oligarchy is 
the point in the evolution of the Solar system where the formation 
of planets transferred from "runaway" growth to "ordered" growth 



4 INTEGRATION AND ANALYSIS 

Expanding on the principle that these "small bodies" only weakly 
gravitationally interact with each other we chose to perform multi- 
ple simulations consisting of the Jovian planets, the Sun and only 
one small body, which is equivalent to a A'-body integration con- 
taining > 2000 objects. Ideally we would have populated the or- 
bital annuli gravitationally stirred by the planets with an oligarch 
number density equivalent to > lOOM®, but these numbers were 
computationally prohibitive. To avoid this complication, we as- 
sumed that oligarch-oligarch interactions can be ignored for N- 
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body simulations run at the Solar system level due to their small 
relative mass. This allowed us to performed a large number of in- 
tegrations containing only one "small body" following which we 
based the remaining analysis on probabilities. At the time of analy- 
sis, this methodology was preferred over a smaller number of sim- 
ulations with a large number of bodies because parallel processing 
code to perform these operations was still being developed. 



4.1 N-Body Integration 

,/V-body integrations were perform using a Burlisch-Stoer inte- 
gration routine e mbedded in the Mercury6 program written by 
IChambersl ( fl99l . The simulation integrates the orbital motions of 
the Jovian planets and one test mass for up to 8 x 10 6 yrs or until 
the test mass either collides with a "big body" (planet or the Sun) 
or is ejected from the Solar system ( > 100 AU from the Sun). The 
accuracy used to define the time step for the integration was de- 
fined to be e = 10~ 15 to ensure accurate results calculated on a 
reasonable time fr ame. Using this pa rameter, integration time steps 
are determined bv lPress et ail d 19971) . 



Hu — H 



l/(2fe+l) 



(1) 



where, H is the user defined time step and tk+i,k is the percent 
difference between the k time steps and k + 1 time steps to span 
H. Mercury6 provides tools to monitor "close encounters" between 
the differ ent objects given the two objects pass within a user defined 
distance dChambersll 19991) . The distance defined to be the minimal 
distance defining a "close encounter" was a half of a Hill radius, 
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Originally the Mercury6 algorithm returned the orbital parameters 
of the two objects that pass within the predefined distance, but we 
altered the code to return the position and velocity vectors of the 
small mass measured with respect to the planet, which were then 
stored for further analysis. 



4.1.1 Initial Conditions 

The orbital parameters used to define the initial orbits of the small 
bodies were all generated randomly with the exception of the ec- 
centricity, which is initially chosen to be equal to zero because this 
is the parameter most sensitive to gravitational perturbation and had 
quickly become randomized during the simulation. In order to re- 
duce the number of simulations that led to stable orbits that do 
not interact with any of the planets, the initial semi-major axis is 
generated by randomly choosing an orbital annulus that lltol J"2002h 
have shown to produce unstable orbits in our Solar system and then 
equated to a random number that would lie within one of these 
rings. Other than the argument of perigee and the ascending node, 
which were chosen to be randomly generated numbers between 
and 2-7T, the angle of inclination was a randomly generated num- 
ber between ±0.01°. Orbital inclinations within this range are well 
within the angular thickness of the Solar nebula after it began to 
"flatten", allow for orbits crossing the planets beyond their orbital 
planes, and are low enough to ensure that the small bodies had or- 
bits crossing the planets. Although planetesimals could have ini- 
tially had large orbital inclinations, objects with inclinations much 
larger than those chosen were less likely to enter a planet's Hill 
sphere. Furthermore, this range of inclinations still allows objects 



to pass through a planet's Hill sphere with essentially a uniform 
distribution as expected. We found (as expected) that the proba- 
bility of locating an object in a "bin" or region around the planet, 
Vbin = (tbin/ttof) has essentially no radial nor angular depen- 
dence. 



4.2 Analysis of iV-Body Integration 

Following the numerical integration performed with Mercury6, a 
more detailed analysis of the encounters between the small body 
and planet(s) was performed. Mercury6 generated vast quantities 
of data consisting of the initial positions and velocities of the small 
body with respect to the planet whenever the small body passed 
within 0.75 Rh of the planet. These data were then used to cal- 
culate the fractional time spent within any given volume element, 
dV, 



dV 



Y[Axi 



around the planet. Due to computational limitations (i.e. RAM as 
well as CPU capacity), dV was restricted to be equal to ~ 10 -7 R% 
corresponding to a Cartesian grid from — 0.75 Rh to 0.75 Rh in 
each direction consisting of 100 bins or divisions setting Aa; = 
0.15 R H . 

The time, dt, spent in each bin for a given pass through the 
Hill sphere is 



ds 

dt= — 

v 



(3) 



where ds is the path of the small body through the volume element 
and both ds and v point in the same direction, dti is approximated 
by assuming that the small body is moving along a straight line 
through the bin and would therefore be equal to 



dti 



(4) 



where Aa^ is the distance to the edge of the volume bin in the 
i th direction where, i — 1, 2, 3 represent the traditional Cartesian 
directions, x, y and z respectively. 

In order to determine the probability of a small body being 
in the j bin, dt is then scaled by the time-of-flight of the small 
body across the Hill sphere, or the amount of time required for the 
small body to pass through our region of interest around the planet 
(~0.75.Rh), 

t tof = J^-(esmhF-F) (5) 

where F is the hyperbolic eccentric anomaly, 

ccshf = e + C ° Sy 
1 + e cos ifi 

and e is the eccentricity, a is the semi-major axis and ip is the an- 
gular measure from periapse. The ratio of dt to ttof provides the 
fractional time spent per encounter in any given volume bin. 

4.2.1 Regularly Spaced Propagation of Trajectory for to 

Determine Where Objects are Likely to be Found within 
the Hill Sphere for Unbound Trajectories 

The trajectory of the small body through a planet's Hill sphere was 
calculated from the angle measured from perigee, <p, and perform- 
ing a coordinate transformation from the orbital coordinates, x', y , 
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to a coordinate system fixed to the planet where the x — y plane is at 
the planets equator and the x axis is directed in the same direction 
as the x axis use for the Solar system. Mathematically, the x — y 
components of the position and velocity vectors in the orbital plane 
are; 



x — r cos if 

y = r sin if 
na 



nil 



■. sin (f 



(e + cos (fi) 



where 



a (1 + e 2 ) 
1 + e cos ip 



(6) 
(7) 
(8) 

(9) 



(10) 



and a is the semi-major axis, e is the eccentricity of the small bodies 
trajectory past the planet calculated from the initial position and 
velocity of the small body as it enters the sphere of interest, and 
/j, = n 2 a 3 . This vector is then rotated to determine where in the 
planet's fixed reference frame the object is located. 

X = n{Q) z pTl(i) x ,n :c sbX' (11) 

where uj is the argument of perigee, Q. is the ascending node and i 
is the angle of inclination and 1Z are the rotation matrices around 
the z axis in the small body's orbital plane, an intermediate x axis 
and the z axis in the planet's orbital plane given by; 




and 



1 

— sin i cos i 
cos i sin i 



TZ(lu) 



cos uj sin u 
— sin ui cos id 
1 



All of the positions and velocities calculated in the orbital frame 
are transformed using equation (II jj , which are then used to deter- 
mine dt and which bin it corresponds to. There were two reasons 
for simulating the orbit in this manner. The most important reason 
was that moving a small body along its trajectory in discrete an- 
gular steps that correspond to the angular width of the bins most 
distant from the planet ensures that the trajectories are sampled 
completely. Since most optimized integration routines will adjust 
the time step to ensure that there is a tolerable offset in a conserved 
quantity like the total energy or angular momentum, it is likely 
that some statistical bins located farthest from the planet would be 
"skipped" over. Therefore the angular step Aip = 1.5 x 10 -2 ra- 
dians was used to correspond to the angular width of a bin at the 
farthest distance of interest, Ax/R ma x, where R ma x = 3/4 Rh- 
To a lesser extent, there is some reduced run-time because Mer- 
cury6 writes data to the hard disk at regular intervals, which would 
then need to be extracted to determine the trajectories. 



5 CALCULATION OF PROBABILITIES 

We calculated the probability of a physical collision between plan- 
etesimals occurring at a given region in space around the planet 



from three specific ratios; (1) the fraction of time spent in the area 
of interest around the planet, (2) the fraction of time spent per bin 
around the planet and (3) the fraction of cross sectional areas (i.e. 
a ratio of the small body cross section to the cross sectional area of 
the region of interest). 

As stated in § 14.2.11 the fractional time spent in the i th volume 
bin is equal to dti /t to f- The fraction of time spent around the planet 
is equal to t p i an ct /tint- For each integration performed, the time a 
small body "spends" around the planet was calculated by totaling 
the time of flight for each "close encounter" for that integration, 



' planet 



(12) 



Mercury6 was configured to return the total integration time, tint, 
for each run since the program would exit if the small body collided 
was ejected of if it collided with a large body (a planet or the Sun). 

The probability that an object is in a given bin, Vbin, provided 
the object is passing through a planet's Hill sphere is, 



T^bin 



dt 

ttof 



planet 



(13) 



Angular dependencies were neglected because Pun varied by less 
than 1 % over 4-zr sr. Although this conveys useful information, it 
is difficult to visualize. More information could be gathered from 
a temporal distribution along a trajectory to determine where in 
the Hill sphere these collisions could occur. Since, dti/t to f has 
little angular dependence, we assumed dti/t to f to be constant 
over all angles, hence an average dti/ttof over 4-7T sr would be 
comparable to averaging multiple temporal distributions to find 
{dt /ttof)r .Therefore, the probability density of finding a small 
body in a given annulus at a given distance (r + dr) per cubic AU 
from the planet, within the planet's Hill sphere is, 



Vr = 



dt 

ttof 



> planet 
tint 



(14) 



The probability that two objects occupy the same bin along a tra- 
jectory through a planet's Hill sphere can be found by integrating 
the probability density V r along an object's trajectory. If these ob- 
jects were large enough to occupy the entire volume of a "bin", 
this would then be the probability of two objects colliding. Unfor- 
tunately this is not the case and this integration must therefore be 
re-scaled by the ratio of the two cross-sectional areas, 



O'sb 
&bin 



Ax 2 



(15) 



If there were only two objects in a planet's Hill sphere at any given 
time, then we would only need the ratio of a a b to <Jbin- Since it is 
likely that there multiple small bodies in the Hill sphere, a s b needs 
to be scaled by the number of small bodies present, A r 'Ppi all et, 
where TV is the total number of small bodies on the annulus that 
interact with the planet and V p i anc t = {tpia.net/tint), leading to a 
cross sectional ratio of 



Csb 
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(16) 



where n(r s b\ Mmmsn) are the number of objects with radius r a b 
and properly scales the area occupied by the "small bodies". The to- 
tal area occupied by the planetesimals will be this ratio, integrated 
over all radii (or at least the radii of interest - 1 to 200 km). Lastly, 
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all of these pieces need to be brought together to yield the proba- 
bility that any two small bodies will collide within the Hill sphere 
of a planet. 



■Pc ol (r) = / VldV^- 

J Chin 

= n(r; MmmsnJ 



(17) 
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Assuming that the rocky material in the Solar nebula is uniformly 
distributed, the surface density of this material for the Solar nebula 
disk between Jupiter and Neptune is 

100 Mm 



OSN 



(18) 



tt(30 2 - 5 2 ) AU 2 

The number density of small bodies orbiting the Sun was calcu- 
lat ed based upon t he siz e dist ribution of planetesimals pre sented 
in INesvornv etall J2003h and ICharnoz & Morbidellil d2003l) . The 
number of small bodies as a function of the MMSN was assumed 
to be 



n(r; A/mmsn) = J)(M M MSN)r" 



(19) 



where 



Mmmsn 
Pe 



r ' V(r)dr 



V(r) is the volume of the planetesimal and p® is the density of 
the Earth. Lastly equation dl7t must be integrated over all "small 
body" radii and the total probability that two planetesimals collide 
within a planet's Hill sphere is 



Vcol = 



Vcoi(r)dr 



(20) 



A/mmsn 



dt 

ttof 



W{f)df 



Ax 2 



-dr 



2 / ^planet 



dV 



where the integration over dV is the volume of the Hill sphere and 
the integration s over both r and f are with r espect to "small body" 
radii g iven by INesvornv et al. ] d2003l) and|ch arnoz & Morbidellil 
d2003h size distribution of dN oc D~ 3 5 dD. 



5.1 Probability of a Bound Orbit and the Resulting Orbital 
Parameters 

The majority of collisions that occurred in a planet's Hill sphere 
did not result in the object remaining gravitationally bound to that 
planet. The magnitude of the kinetic energy of a typical collision is 
extremely high (figure[T} and is large enough to completely destroy 
objects with radii < 20 km. Furthermore, even if the two objects 
were able to coalesce, the majority of collisions did not remove 
sufficient orbital energy from the combined object to result in a 
bound orbit. 

However, to simplify matters, we can make some initial as- 
sumptions about the dynamics of a collision needed to be made in 
order to determine which collisions would result in a bound orbit 
or not and revisit these assumptions given the results. The first as- 
sumption was that the two objects would "stick" together, or 



vi — v 2 < \1Gpr 



(21) 



where the objects have a density, p s j,, equal to the density of the 
Earth, p® w 5gcm -3 . We understand that this is not the case, 



but calculating the orbital characteristics of the resulting inelastic 
collision is simpler and can serve as a lower limit for the proba- 
bility of objects being captured. Although technically not an as- 
sumption, the trajectories known to pass through a given volume 
bin (d < 0.01 AU, where d is the distance between the two objects) 
were "forced" to collide by translating the small bodies to the av- 
erage position vector with their given velocity. This translation re- 
sults in a error of < 1%, and the probability of a collision occurring 
was already determined in equation [20] Furthermore the ensemble 
required to produce even one collision is prohibitively large, "forc- 
ing" trajectories to "cross paths" is a reasonable approximation (as- 
sumption). Lastly, provided equation d2U is satisfied, a collision 
results in a bound orbit provided, 

. GM planet 



Va-m < 2- 



Tcol 



(22) 



where V cm = (mi vi + m^vi)! (mi + U12) is the center of mass 
velocity, r co i = (ri + ra)/2 is the point where the collision is 
"forced" to occur and that the agglomerate doesn't smash into the 
planet (i.e. a(l - e) > r planet). 

Initial conditions for trajectories corresponding with all of the 
"close encounters" for each planet were stored so that we could use 
each unique initial condition to integrate trajectories and determine 
collision characteristics. Limiting ourselves to 500 unique initial 
conditions, we assign an initial condition to one object, then step 
through the remaining conditions counting the total number of "as- 
sumed" collisions (d < \^2Ax), N c , and the number that result 
in a bound orbit, iVj, (resulting object satisfying equationl22l. This 
results in 



(i - 1) ~ 4 x 10 7 



integrations. These integrations produced ~ 10 5 collisions of 
which < 10 resulted in a bound orbit (table ??). 

Figure ((2} are plots of the eccentricity as a function of semi- 
major axis for an irregular satellite captured by the planets via an 
inelastic collision between two planetesimal. This plot shows that 
the eccentricities of orbits arising from these collisions are essen- 
tially random, although slightly skewed with a mean of ~ 0.7. The 
region above the line corresponds to orbits with periapse less than 
or equal to the radius of the planet (in this case Jupiter). These or- 
bits would result in a collision with the planet. Additionally, figure 
10 shows the distribution of inclinations from the resulting bound 
orbits. These distributions are essentially random distributions with 
a couple of the planets having a slightly higher probability of orbits 
with marginally retrograde orbits with inclinations ~ 120 deg. 



6 RESULTS 

We ran 2000 random simulations resulting in > 2 x 10 8 close en- 
counters between a given planet and a test mass that revealed the 
mean fractional times, (t p i a net /tint) spent around Jupiter, Saturn, 
Uranus and Neptune are 5.48 x 10~ 5 , 5.32 x 10" 5 , 1.31 x 10~ 5 and 
6.29 x 10~ J respectively (described in §0. As stated in § 14.11 nu- 
merical integrations simulated periods of ~ 8 Myrs unless the test 
mass was ejected or collided with either a planet or the Sun prior 
to the 1 Myr period being completed. Although we attempted to 
place the test mass es o n known unstable orbits around the Sun de- 
termined from Ito (2002j), ~ 40% of the simulations ran the entire 
8 Myrs, ~ 25% ended < 1 Myr and the remaining 35% uniformly 
distributed between 1 — 5 Myrs. 
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Figure 1. Ratio of collision energy j5 C o1 to gravitational binding energy Bbind as a function of the "small body" (planetesimal) radius. The (red) dashed lines 
represent one a error for the measured collision energies. The (red) dashed lines correspond to the measured uncertainty of E co \ 



Table 1. Collision Rate of Irregular Satellites (Myr 1 ) 



Object Radius (km) 


Jupiter 


Saturn 


Uranus 


Neptune 


1 < r < 10 


0.601 ±0.304 


40.30 ± 1.25 


133.1 ±4.2 


176.3 ± 11.0 


10 < r < 100 


(1.90 ±0.11) x 10" 3 


(1.28 ± 0.04) X 10" 1 


(4.21 ±0.13) X 10" 1 


(5.57 ±0.35) x 10" 1 


100 < r < 1000 


(6.01 ±0.34) x 10" 7 


(4.03 ±0.01) x 10" 3 


(1.33 ±0.04) x 10" 3 


(1.76 ±0.11) x 10" 3 



Using the formulae presented in § [5] we found that the frac- 
tional time that a test mass spends around any given planet was 
~ 10~ 5 , therefore if an object remained gravitationally bound to 
the Sun on an unstable orbit for 1 Myr it would be likely to spend 
~ 50 yr within the Hill sphere of any one of Jovian planets. As- 
suming that the MMSN is ~ 100 M®, collision rates for varying 
object sizes are presented in table ([T}. For object radii between 1 
and 10 km, collision rates are 0.601 ± 0.304 Myr _1 for Jupiter, 
40.3 ±1.25 Myr _1 for Saturn 133.1 ±4.2 Myr" 1 for Uranus and 
176.3 ± 11.0 Myr" 1 (figure|3j for Neptune. Although these col- 



lision rates may appear to be large and imply that there should be 
a significantly larger population of irregular satellites, when com- 
bined with the probability that a collision between two planetes- 
imal results in a bound orbit, capture rates for objects with radii 
between 1 and 10 km are (table |2) (2.00 ± 0.22) x 10~ 5 Myr" 1 , 
(6.74 ± 0.64) x 10" 4 Myr" 1 , (4.19 ± 0.46) x 10" 3 Myr" 1 , 
(6.40±0.52) x 10" 3 Myr" 1 , for Jupiter, Saturn, Uranus and Nep- 
tune respectively, for a MMSN of ~ 100 M®. Given energies plot- 
ted in figure [T] it is possible that a collision between an object with 
radius ~ 10 km collided with one with a radius of ~ 100 km, 
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Figure 2. Eccentricity plotted as a function of Semi-Major axis for the conglomerate generated by the collision of two objects within Jupiter's Hill sphere. 
The region above the (red) line correspond to the bound orbits that pass within Jupiter's radius thereby resulting in a collision with the planet and the vertical 
(green) line the semi-major axis of the "regular" moon orbiting farthest from the planet. There are no restrictions placed on the collision energy for the points 
plotted above and the average value for V 2 , is ~ 4 X 10 3 H a b/ r coh the binding energy per unit object of the small body where we assumed r co i ~ 10 km. 
The different points correspond to different object ratios between the colliding objects: (black) triangle are 1 to 1, (red) diamonds are 10 to 1, (green) pentagons 
are 100 to 1 and (blue) hexagons are 1000 to 1. 



which would have resulted in the larger of the two being broken 
apart. Another point of interest is that Jupiter has a "capture" rate 
significantly lower than the other planets. This is most likely due 
to the fact that as these objects move closer to the Sun, their speed 
increases, thereby requiring more energy to be removed when the 
two planetesimals collide in order to capture the satellite. 



7 CONCLUSIONS 

The Jovian planets are known to have between ~ 20 and ~ 60 
irregular satellites orbiting them depending on the planet (Jupiter 
> 60, Neptune ~ 10). Assuming that these satellites were all cap- 
tured over a period of 1-2 Myrs, the "capture rates" for the plan- 
ets would be ~ 0.001 Myr . In § [6] we present the probabil- 



ity that a collision resulted in a bound orbit Pbound ~ 10 -5 per 
Myr for a total MMSN containing about 100 and object radii 
1 km < r < km. The collision energy scaled by the binding en- 
ergy of an object with a 100 km radius (figure [T} shows that col- 
lisions between the different objects would have to have a mass 
ratio < 0.1 to ensure that both weren't completely destroyed and 
a ratio of ~ 1 to be broken apart and not reduced to rubble. 
Using a MMSN based upon asteroid measurements presented by 
iMorbidelli et alj d2009|). which is timmsn ~ 10 6 AU~ 2 and size 
distri bution given in dNesvornv et al. I l2003l : Ich arnoz & Morbidellil 
120031) . dN oc D 3 5 dD, the rate the Jovian planets could expect 
to capture irregular satellites through collisions is ~ 10~ 3 Myr~ 



(figure |3). This would be consistent with the iGoldreich et al.l 
(2004b) oligarchy models with the exception that objects primarily 
consist of > 10 — 100 km sized objects, not ~ 1 km and requiring 
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Figure 3. Collision rates for small body collisions occurring within the planet's Hill sphere. The number of small bodies from calculated using equation 119) 
is used in conjunction with equation (20} to calculate collisions per Myr. The (red) dashed lines correspond with the propagated measurement uncertainty of 

{dt/t to f) r and (t p i anct / 'tint) 



Table 2. Capture Rate of Irregular Satellites (Myr -1 ) 



Object Radius (km) 


Jupiter 


Saturn 


Uranus 


Neptune 


1 < r < 10 
10 < r < 100 
100 < r < 1000 


(2.00 ±0.10) X 10" 5 
(6.32 ±0.69) x 10~ 8 
(2.00 ±0.22) x 10 -10 


(6.74 ±0.64) x 10" 4 
(2.13 ±0.20) x 10~ 6 
(6.74 ±0.64) x lO -9 


(6.74 ±0.64) x 10" 4 
(1.33 ±0.15) x 10~ 5 
(4.19 ±0.46) x 10~ 8 


(4.19 ±0.46) x 10~ 3 
(2.02 ±0.17) x 10~ 5 
(6.40 ±0.52) x lO -8 



a slightly larger MMSN in order to populate the Edgeworth-Kuiper 
belt, asteroid belt and Oort cloud. Since the probability is highly 
dependent upon the size distribution of planetesimals as well as the 
mass of the protoplanetary disk, we also invest igated size distribu- 
tions of dN oc D 18 dD jgottke et al ]l2010bl) . dN oc D~ 2 2 dD 
and dN oc D~ 2 S dD dTerai & Itoh]|2010h all producing similar 
results for object radii 10 km < r < 100 km (e.g. same order 
of magnitude). However for 1km < r < 10 km, collision rates 



for a size distribution dN oc D l s dD is approximately half that 
for dN oc D~ 3 - 5 dD and for the radii 100 km < r < 1000 km 
its approximately doubled. As the magnitude of the power law in- 
creases, naturally the probability of collisions increases for smaller 
objects because there are a larger number of them. Furthermore, 
allowing plane t s to m igrate as in the Nice model or simulated by 
iLvkawka et all (2009), each planet would sweep through a larger 
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Figure 4. Distribution of the angle of inclination for the bound orbits represented in figure [2] As stated in the assumptions, these inclinations are the result of 
a collision between to planetesimals within a Jovian planets' Hill sphere where the two objects "stuck" together. 



effective area, thereby marginally increasing the probability of a 
collision between objects. 

The results presented above are subject to some caveats. 

• The two objects colliding, "stick" together and form a single 
object as a result of the collision, which can not be the case because 
the collision energies are > 10 times the gravitational binding en- 
ergy. This assumption can result in an underestimation in the orbital 
energy removed from the collision and lead to an underestimation 
of the number of bound satellites. 

• The objects are assumed to never achieve a "quasi-bound" or- 
bit, which is known to occur where objects can remain "bound" to 
the planet for more than 100 years. Although these objects are rare, 



the amount of time spent within the Hill sphere is > 10 4 times those 
scattering off of the planet. Furthermore since these objects are al- 
ready on a "quasi-bound" orbit around the planet, less energy needs 
to be removed from the collision. Each of the impacts from this as- 
sumption, like the first, would lead to an underestimat e of the num- 
ber of bound satellites. iHorner & Wvn Evansl d2006l) have shown 
that the Centaurs can be capture like Trojan satellites and even ir- 
regular satellites are able to be temporarily captured by Jupiter for 
> 10 4 years. Also, in recent yea rs a number o f come ts, Gehrels, 
Shoemaker-Levy 9 and Oterma Ohtsuka et al. (2008) and refer- 
ences therein. Furthermore. IHorner & Lvkaw ka (2010) have shown 
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such captures to be quite frequent. Such events could significantly 
boost our capture rates. 

• All of the planets formed at there present location. It has al- 
ready been accepte d that, at least, Neptune has migrated outward 
since its formation jMalhotrall 993). As planet migrated away from 
the Sun, there Hill sphere increases, which would mean that for 
planets moving outward (Neptune and Uranus), their Hill spheres 
would have been over estimated, hence the probability of capture 
as well. We believe that this over estimation is compensated by the 
fact that as the planet(s) migrated the annulus swept out by their 
orbit(s) grows giving them "access" to a larger number of objects. 
However, for planets that migrated inward, their Hill spheres were 
underestimated, this coupled with the "growth" of their orbital an- 
nulii would lead to an underestimation of the number of bound 
satellites. 

• The density of objects colliding is ~ p®, where in actuality 
the density is less than the density of the Earth because the Earth 
and Mercury are the objects with the greatest known density. This 
will lead to an underestimate in the number of objects (planetesi- 
mals), hence produce an artificially low collision rate. 



We found that collisions between objects with a mass ratio as large 
as 1000 to 1 only marginally reduce the capture rate (probability of 
bound orbits listed in table ??). It is possible for a large mass to col- 
lide with a smaller mass and, though the smaller of the two would 
be destroyed because the collision is large enough to reduce the 
small mass to rubble, the larger mass would fare much better and 
would most likely simply be broken apart, resulting in some of the 
debris from the collision remaining bound to the planet. Further- 
more, since ~ 1000 1 km radius objects are contained in a mass 
with a 100 km radius, it is possible that only a small fraction of 
the remnants from the collision could remain bound to the planet 
were the remainder did not have enough energy dissipated by the 
collision. A detailed analysis of the collisions of these objects and 
simulations of the distribution would need to be done to determine 
what fraction of these remnants would remain bound to the planets. 

We found that the average semi-major axes for bound orbits 
are essentially randomly distributed between 0.3 and 2 Hill radii, 
which is consistent with observation with the exception that many 
of the objects captured in the simulation had highly eccentric orbits 
and many with a semi-major axis > 5 Hill radii, which is not what 
has been observed. This is most likely because these objects were 
probably stripped from the planet by the Sun's gravity or interacted 
with the regular satellites to produce more circular orbits, closer 
to the planet. The act of satellites being stripped from the planets' 
orbit is yet another factor not accounted for in the analysis above. 
Since nearly half the capture objects have orbits that are on the 
outer cusp of the Hill sphere, either their orbits would have to decay, 
or they would be stripped from the planet. 

Although it is unlikely that all of the irregular satellites were 
captured via collisions with the host planet's Hill sphere, we have 
shown that there is a finite probability that irregular satellites can 
be capture through collisions, which would serve as a lower limit 
based on the assumptions that have been made. Additionally, the 
large collision ener gies that we measured w ould be consistent with 
the observations of INesvornv et alj d2004l) who state that some of 
Jupiter's irregular satellites can be clustered into specific groups 
with similar orbital elements and spectral characteristics. 



8 FURTHER WORK 

There are several issues that need to be addressed in future work, 
which can be subdivided into two different sections. One section 
incorporates the assumptions made in the simulations performed 
and the other the short comings of the current model. 

The condition s modeled in thi s resear ch most closely mimic 
those presented by Goldrei ch et al.l (|2004a) where the planets are 
formed at their current orbital locations. It differs in the sense that 
the energy of planetesimal orbits is not dissipated, "cooled", and 
steadily increases throu gh interactions with the planets, "heating". 
Goldre ich et al. (2004a) state that interactions between planetesi- 
mal will effectively remove much of this energy by dispersing it 
amounts the entire planetesimal population referred to as "cooling". 
Planetesimal cooling will effectively reduce the energy of these ob- 
jects and the resulting collision energy, therefore should increase 
the probability of a bound state. 

In addition to accounting for orbital " cooling", dependence 
on t he migration of p lanet orbits presented in lTsigapis et alj( 2005) 
and Lykawka et al. I <2009h must be investigated. INesvornv et al.l 
(2007) present a mechanism of capture based upon a three body 
interaction between to gas giants and a planetesimal passing be- 
tween them during a close encounter. This theory can recreate the 
observed orbits for irregular satellites around Saturn, Uranus and 
Neptune, but fails for Jupiter since th ere are no known case s of 
Jupiter encountering the other planets. INesvornv et al.l d2007l) go 
on to state that Jupiter would have to capture its satellites through 
some other means. It would be beneficial to determine if this mi- 
gration increases collision rates by increasing the effective area of 
the feeding zones for the planets, hence the number objects passing 
through the planets' Hill s pheres as well as prod uce bound orbits 
similar to those observed bv lNesvornv et al.| (2007) and comparable 
to observed orbits. 

Lastly, due to the high magnitude of the collision energy and 
low probability of a collision between two objects resulting in a 
bound orbit, it is important study the dynamics of the collisions 
themselves. Since it is most likely that the objects are broken apart 
during the collision, the kinematics of these inelastic collisions and 
the kinetic energy of the remnants should be carefully studied in as- 
teroid families to determine the probability that enough energy has 
been removed from the remnants of the collision to remain bound 
to the planet. It is likely that this would result in a much higher 
probability of bound orbit. 
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